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ABSTRACT 

We present almost-simultaneous detections of Cygnus X-1 in the radio and mm 
regimes, obtained during the low/hard X-ray state. The source displays a flat spectrum 
between 2 and 220 GHz, with a spectral index |a| < 0.15(3tT). There is no evidence 
for either a low- or high-frequency cut-off, but in the mid-infrared (~ 30/im) thermal 
emission from the OB-type companion star becomes dominant. The integrated lumi- 
nosity of this flat-spectrum emission in quiescence is > 2 x 10'^^ erg s~^ (2 x 10^'' W). 
Assuming the emission originates in a jet for which non-radiative (e.g. adiabatic ex- 
pansion) losses dominate, this is a very conservative lower limit on the power required 
to maintain the jet. A comparison with Cyg X-3 and GRS 1915-1-105, the other X-ray 
binaries for which a flat spectrum at shorter than cm wavelengths has been observed, 
shows that the jet in Cyg X-1 is significantly less luminous and less variable, and is 
probably our best example to date of a continuous, steady, outflow from an X-ray bi- 
nary. The emissive mechanism reponsible for such a flat spectral component remains 
uncertain. Specifically, we note that the radio-mm spectra observed from these X-ray 
binaries are much flatter than those of the 'flat-spectrum' AGN, and that existing 
models of synchrotron emission from partially self-absorbed radio cores, which predict 
a high-frequency cut-off in the mm regime, are not directly applicable. 
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1 INTRODUCTION 

Cygnus X-1 (V1357 Cygni, HDE 226868) is one of the bright- 
est X-ray binaries and the classical black hole candidate 
(BHC); the inferred mass for the compact object is ~ 7 
M0 (e.g. Gies & Bolton 1986). The radio emission in the 
low/hard X-ray state in the cm-wave band is persistent, 
and although it varies by a factor of about 5 the varia- 
tions are much less spectacular than those of other well- 
studied XRBs such as Cygnus X-3 and GRS 1915-f 105. The 
radio flux is modulated at the 5.6-day orbital period, par- 
ticularly at higher frequencies, and has a spectral index 
(a = Alog S^/Alog i^) close to zero (|Qf| < 0.1 between 2 
- 15 GHz; Pooley, Fender & Brocksopp 1999). There are 
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also variations on timescales as short as 1 hour and long- 
period fluctuations with an apparent 140-day period. The ra- 
dio emission is also known to change at major state-changes 
in the X-ray emission (Hjellming, Gibson & Owen, 1975), in 
the sense that it is reduced in the X-ray high/soft and/or 
'intermediate' states. For long periods both the radio and 
X-ray emission are relatively stable, with the mean radio 
flux density at cm wavelengths in 1996 - 1998 being about 
14 mjy (Pooley et al. 1999). A detailed comparison of short 
and long-term behaviour from radio to hard X-rays is pre- 
sented in Brocksopp et al. (1999). Stirling, Spencer & Gar- 
rett (1998) and de la Force et al. (in prep) have resolved 
the radio emission from Cyg X-1 on milliarcsecond angular 
scales, supporting an origin in a spatially extended outflow 
or jet. A single previous detection of 10 ±3 mJy at 250 GHz 
reported by Altenhofi', Thum & Wendker (1994), which was 
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Table 1. Log of IRAM observations of Cyg X-1 in 1997 Aug 



Table 2. Log of JCMT observations of Cyg X-1 in 1998 May 



Date Frequency Flux density 
(MJD) (GHz) (mjy) 



50644 



15.9 ±4.9 



not simultaneous with radio observations, implied the pres- 
ence of an approximately flat spectrum from cm through 
mm wavelengths. 

In this paper we describe multiple detections of the 
source simultaneously at cm and mm wavelengths. These 
observations reveal a very flat spectrum extending from the 
radio to the mm regimes. They do not reveal any high- 
frequency cut-off in the spectrum and raise the inferred qui- 
escent cm-mm luminosity by a factor of about 15. 



2 OBSERVATIONS 

2.1 Radio 

Cygnus X-1 is monitored approximately daily by the Green 
Bank Inteferometer (GBI) at 2.3 & 8.3 GHz and by the 
Ryle Telescope (RT) at 15 GHz. Details of typical GBI ob- 
servations can be found in e.g. Waltman et al. (1994). The 
RT observations of Cyg X-1 are described in Pooley et al. 
(1999). 



2.2 IRAM 

The observations were carried out with the IRAM 30m tele- 
scope at Pico Veleta (Spain) on 1997 Aug 4. We observed 
Cyg X-1 simultaneously in the HCN (J=l-0) transition at 88 
GHz, with the SIS 3mm-l receiver, in CN (N=l-0) transition 
at 113 GHz, with the SIS 3mm-2 receiver, in CN (N=2-l) 
transition at 220 GHz using the R.230G1 (1mm) receiver. 
The half power beam widths (HPBW) are respectively 27" , 
22", 10.5" at 88 GHz, 113 GHz and 226 GHz. The observa- 
tions were taken using a position switching mode, with a ref- 
erence position at -1 arc min in RA compared to the source 
position. The mean system temperatures were 185K, 403K 
and 514K, and the opacities 0.1, 0.28 and 0.51 respectively 
for the 3mm-l, 3mm-2 and 230G1 receivers. Due to the fac- 
tor 3 between the opacity at 88 Ghz and 113 GHz, we used 
the 3mm-l receiver, with a 512 MHz bandwidth, centered on 
the HCN (J=l-0) line to estimate the Cyg X-1 continuum 
emission. The calibration was taken on the source K3-50A, 
and we observed Cyg X-1 for a real time of 30 min, using 
12 subscans of 20 s (10 s on source and 10 s off source). The 
continuum emission from Cyg X-1 was detected at a level of 
15.9 mJy (3.2a) assuming a 6 Jy/K ratio at 88.6 GHz 



2.3 JCMT 

JCMT SCUBA (Holland et al., 1999) 1350 and 2000 ^m pho- 
tometry observations of Cyg X-1 were carried out in 1998 
May. Rather than employing the arrays these observations 
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146 
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all data 


146 


8.9 ±2.1 


- 


221 


13.3 ±1.6 



- t 

1997 July 24 


t t * 

- Aug 15 (MJD 50653 - 50675) 


. 


t 
1998 May 10 - 


t t 

- June 3 (MJD 50944-50966) 


* 

t 1 



Frequency (GHz) 

Figure 1. Mean 2 — 220 GHz radio/mm spectra of Cyg X-1 in 
1997 August and 1998 May. The spectrum is clearly flat across 2 
decades of frequency at both epochs (best-fit spectral indices of 
0.07 ± 0.04 and -0.06 ± 0.05 in 1997 and 1998 respectively). 



used the longer wavelength single-pixel bolometers posi- 
tioned around the arrays. The data reduction was performed 
in the standard manner using SURF (Jenness, 1998; Jenness 
& Lightfoot, 1998). Photometric calibration was achieved 
by skydip analysis and photometry of Uranus and/or the 
secondary standard CRL2688. The emphasis of the obser- 
vations was on 'detection' rather than accurate flux deter- 
minations; the 1350 /jm flux densities have an accuracy of 
about 20%, the 2000 ^m flux densities of 25-30%. The radio, 
IRAM and JCMT data are plotted in Fig 1. Least-squares 
fitting of single power laws to the data from the two epochs 
results in spectral indices of 0.07 ±0.04 and —0.06 ±0.05 for 
1997 August and 1998 May respectively. 



2.4 XTE 

Cyg X-1 is monitored up to several times daily in the 2-12 
keV band by the Rossi X-ray Timing Explorer (RXTE) All- 
Sky Monitor (ASM). See e.g. Levine et al. (1996) for more 
details. The total flux measured by individual scans is plot- 
ted in the top panels of Figs. 2. Cyg X-1 was in the low/hard 
X-ray state throughout the period of these observations. 



3 DISCUSSION 
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Table 3. Mean 2 - 220 GHz flux densities of Cyg X-1 in 1997 
August and 1998 May. 'Errors' at t" > 89 GHz are dominated 
by measurement uncertainties, whereas those at 2-15 GHz are 
dominated by intrinsic source variabiUty. These data are plotted 
in Fig 1. 



Cygnus X-1 1997 Jul 25 - 1997 Aug 14 



Date 
(MJD) 



Frequency 
(GHz) 



Flux density 
(mjy) 



50664 - 


89 


15.9 ±4.9 


50675 


15 


11.5 ±4.0 


(1997 July - 


8.3 


12.8 ±4.6 


August) 


2.3 


12.5 ±4.0 


50945 - 


221 


13.3 ±1.6 


50964 


146 


8.9 ±2.1 


(1998 May - 


15 


10.5 ±3.0 


June) 


8.3 


14.0 ±6.0 




2.3 


14.4 ±4.2 



3.1 Flux densities at the two epochs 

Table 3 summarises the mean flux densities observed dur- 
ing these two epochs. At mm wavelengths (frequencies > 89 
GHz) errors reflect measurement uncertainties; at cm wave- 
lengths 'errors' reflect the intrinsic source variability, about 
which little is known at mm wavelengths. 

3.2 Variability at cm and X-ray wavelengths 

The variations of 15-GHz flux density and X-ray count-rate 
during the two observing intervals concerned are shown in 
Fig 2. Both intervals show some enhanced activity in the 
X-ray band; see discussion in Brocksopp et al. (1999). The 
overall mean amplitude of the 5.6-day modulation at 15 GHz 
is 2.0 mJy (zero-peak). Table 3 lists the r.m.s. variability 
at 2, 8 & 15 GHz over a period of 22 days (i.e. ~ 4 orbital 
periods) centred on each of the two mm observation periods. 
In the light of these variations, and the absence of exactly 
simultaneous observations, it is possible that the spectrum 
in the cm - mm regime also has some short-term variations, 
but it is probable that the mean spectral index is close to 
zero up to 220 GHz. 



3.3 The broad band radio — optical spectrum and 
a comparison with other X-ray binaries 

In Fig. 3 we show an extrapolation of the flat radio - mm 
spectrum through the infrared spectral region, combined 
with published infrared flux density measurements. It is 
clear that shortwards of ~ 30/im, the observed emission will 
be dominated by the thermal component from the OB-type 
companion star and its wind. Clearly in Cyg X-1 even if the 
flat spectrum emission did extend to ~ l/im it would be 
extremely hard, probably impossible, to detect, being more 
than two orders of magnitude weaker than the thermal emis- 
sion. 

Only two other X-ray binaries have been detected with 
a flat spectrum extending from cm to shorter wavelengths. 
In GRS 1915+105 the flat spectrum (inferred synchrotron) 
oscillations are observed from 13 cm to 2/im (e.g. Fender et 
al. 1997; Fender & Pooley 1998). In Cyg X-3 the flat spec- 
trum is observed to 0.85 mm (Fonder et al. 1995; Ogley et al. 
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MJD-50000.0 
Cygnus X-1 1998 May 10 - 1998 Jun 01 



Figure 2. The variations of the 2-12 keV X-ray count rate (upper 
panels) and the flux density at 15 GHz (lower panels) over the 
two intervals during which mm-wavc observation were made. 



in prep), and maybe also to 2/im (Fender et al. 1996). Both 
of these systems are generally brighter radio sources than 
Cyg X-1, and yet also more distant. Table 4 compares the 
cm-mm luminosities of the three sources for periods when a 
flat spectrum is detected (i.e. Cyg X-1 in the low/hard X- 
ray state, GRS 1915+105 during periods of oscillations, Cyg 
X-3 almost all the time). It is clear that the flat-spectrum 
component in Cyg X-1, as well as being less variable, is also 
considerably less luminous than those observed from Cyg 
X-3 and GRS 1915+105, by at least an order of magnitude. 
As noted already in Fender et al. (1997) the flat-spectrum 
emission in Cyg X-3 and GRS 1915+105 appears to have 
approximately the same luminosity. 

GX 339-4 is a persistent black hole candidate X-ray bi- 
nary with similar radio properties to Cyg X-1 (Hannikainen 
et al. 1998; Fender et al. 1999 and references therein). In par- 
ticular the source displays at cm wavelengths a flat spectrum 
with comparable luminosity to that of Cyg X-1. We fully ex- 
pect therefore that sufficiently sensitive observations should 
also detect a flat spectrum through the (sub)mm regime 
from this source. Additionally, as GX 339-4 is believed to 
be a low mass X-ray binary with a less luminous companion 
star than in the Cyg X-1 system, we may have more chance 
of detecting the flat spectrum at near-infrared wavelengths. 
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. ISOCAM data (Mirabel et al. 1996) 

* IR Phometry (Persi et al. 1980) 




Prequency (Hz) 

Figure 3. The broadband spectrum of Cyg X-1 from radio to 
optical wavelengttis, combining our radio/mm data with the ISO- 
CAM observations of Mirabel et al (1996) and ground-based pho- 
tometry of Persi et al. (1980). It is clear that even without a high- 
frequency cut-off the flat spectral component will be dominated 
by the thermal emission from the OB supergiant at wavelengths 
of 30/im or shorter. 



Table 4. Comparison of cm— mm luminosities of Cyg X-1, Cyg 
X-3 and GRS 1915-1-105 

Source 5,y(flat) assumed relative 

(mjy) distance (kpc) luminosity 



Cyg X-1 

Cyg X-3 

GRS 1915+105 



10-15 
50-100 
30-60 



2-3 
7-10 
7-12 



20-250 
30-215 



3.4 The nature of the flat-spectrum component 

Tables 1-3 and Figs 1 & 3 summarise the observations of 
Cyg X-1 for 1997 Aug 4 and 1998 May 11-20. The source 
is clearly displaying a flat spectrum through the radio-mm 
regimes at both epochs. While radio emission from X-ray 
binaries is generally assumed to be synchrotron in origin 
(see e.g. Hjellming 1988; Hjellming & Han 1995), in the case 
of Cygnus X-1 we do not have direct observational evidence 
for this. Even the most rapid variability observed at 15 GHz 
does not require a brightness temperature in excess of 10® 
K, and there is no direct measurement of linear polarisation. 
So, while some form of self-absorbed synchrotron emission 
remains a possible origin for the flat spectral component, 
other emissive mechanisms must also be considered. 



3.4-1 Energetics 

The observed luminosity of a flat-spectrum source is directly 
proportional to the total bandwidth. In the case of Cyg X-1, 
the cm-mm flat spectral component corresponds to a radia- 
tive luminosity of > 2 x 10^^ erg s"^ (2 x 10^'' W). If the 
emission arises in an outflow in which non-radiative (e.g. 
adiabatic expansion) losses dominate (which seems likely to 
be the case for relativistic jets from X-ray transients, see 
e.g. Hjellming & Han 1995) then even the integrated ra- 
diative luminosity is only a lower limit on the total power 
(i.e. it neglects e.g. electron acceleration and bulk kinetic en- 



ergy) required to maintain the jet. Beyond the mm regime, 
in the infrared, thermal emission from the companion, stel- 
lar wind and accretion disc begin to dominate the spectrum 
of the system (Fig. 3) and it may be very difficult to ever 
measure any high-frequency limit to the flat spectral com- 
ponent emission. Note that there is strong observational 
evidence that the flat-spectrum oscillations observed from 
GRS 1915+105 are dominated by adiabatic expansion losses, 
based upon the similarity of the oscillation decay rates at 
cm and infrared wavelengths (Fender et al. 1997; Fender & 
Pooley 1998). 



3.4-2 Partially self-absorbed synchrotron models developed 
for 'flat-spectrum' AGN 

It is easy to draw parallels between the flat-spectrum radio- 
mm emission from Cyg X-1 (and also Cyg X-3 and GRS 
1915+105; see above) and the 'flat-spectrum' extragalactic 
radio sources. These systems are generally radio-loud AGN 
in which the flat-spectrum component corresponds to the 
'core' or base of the jet. As pointed out by Cotton et al. 
(1980) it would appear to require a 'cosmic conspiracy' of 
superposition of individual self-absorbed synchrotron com- 
ponents in order to produce a composite flat spectrum. 
Marscher & Gear (1985) and O'Dell et al. (1988) showed 
that you can more comfortably reproduce 'flat-spectrum' 
variability via shocks in conical jets. A conical jet model for 
radio emission from X-ray binaries was presented by Hjellm- 
ing & Johnston (1988). Giovanoni & Kazanas (1990) sug- 
gested that energy transport by relativistic neutrons natu- 
rally explained the combination of electron spectrum, den- 
sity and magnetic fleld profiles required to produce an ob- 
served flat synchrotron spectrum. Alternatively, Wang et 
al. (1997) have suggested that the flat-spectrum emission 
is optically thin from a flattened electron energy distribu- 
tion. However, there are problems with the application of 
most, possibly all, of these models to the flat radio-mm(- 
infrared) spectra observed from Cyg X-1, Cyg X-3 and GRS 
1915+105. The simultaneous radio-infrared oscillations ob- 
served in GRS 1915+105 constitute evidence against both 
shocks which cool via radiative losses (as the decay rate is 
the same at 2 cm and 2 /xm) and an optically thin solution 
(as the infrared-radio delay, as well as delays within the 
radio band, suggest significant optical depth effects). Fur- 
thermore, all the conical jet and related models only predict 
a flat spectrum over at most three decades in frequency; the 
problem in all cases is the prediction of a high-frequency cut- 
off somewhere in the mm band. This is observed in nearly 
all cases for 'flat-spectrum' AGN, where the mean spectral 
index in the mm band is in fact < Omm >= —0.75 ± 0.05 
(Bloom et al. 1994). It is therefore clear that the three X- 
ray binaries in question have much flatter (consistent with 
completely flat) radio-mm(-infrared) spectra than the 'flat- 
spectrum' AGN, and the applicability of the self-absorbed 
synchrotron models to these X-ray binary spectra remains 
to be established. 

If the emissive mechanism is synchrotron, then as- 
suming that the mm emission is not signiflcantly Doppler 
boosted, we can estimate a minimum size for the emitting re- 
gion from the inverse Compton brightness temperature limit 
of lO^^K. At 220 GHz, this is only 10^" cm, which is rela- 
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tively close to the compact object and well within the binary 
separation (~ 10^^ cm) of the system. 

3.4-3 Alternatives to synchrotron emission? 

An obvious candidate for the emissive mechanism of a flat 
spectral component is optically thin free-free emission. For 
optically thin free-free emission from a thermal plasma, we 
need to have a sufficiently large emission measure whilst 
keeping the spectrum optically thin to u <2 GHz. Assuming 
a fully ionised pure hydrogen plasma and a Gaunt factor of 
unity (neither of which assumptions will affect an order-of- 
magnitude estimate), and a distance to the system of 2.5 
kpc, we find that we need to satisfy the following criteria: 

r'^NlT^^''^ > 4 X 10^** 

and 

rNlT~'^''^ < 2 X 10^" 

where r is the dimension of the cloud (cm) along the line 
of sight, A'^e is the electron number density (cm~^) and Te 
is the electron temperature (T). The first criterion is neces- 
sary to produce the observed level of emission, the second to 
prevent the cloud becoming optically thick to free-free self- 
absorption. As a result we can determine a minimum size 
of a (spherical) cloud (and corresponding Ne) for different 
temperatures. For a cloud of T = 10*K, i.e. in approximate 
thermal equilibrium with the OB star wind, r > 10^** cm 
(A^e ~ 10® cm~^). For a much hotter cloud of temperature 
lO'^K a dimension of r > lO^"* cm (A^e ~ 10^ cm"^) is still 
required. This is very large indeed compared to the dimen- 
sions of the binary orbit, and a significantly larger emission 
measure than would be expected for the OB star alone. In 
this case the 50% variability timescale, would be > 1 hr 
for a lO^K cloud, and > days for T = 10*K. The small 
(1-1-4.4 X 10~^*'T) correction for relativistic free- free emission 
is insufficient to significantly alter the result. Nonthermal 
optically thin free-free emission should also produce a fiat 
spectral component, with (potentially) a greater emissivity 
than thermal free-free emission, but precise determination of 
this (including calculation of the relevant nonthermal Gaunt 
factors) is beyond the scope of this paper. Regardless, as 
noted above it is difficult to invoke a purely optically-thin 
solution as there is evidence for frequency-dependent delays, 
indicating a significant optical depth. 

Wright & Barlow (1975) have calculated the spectrum 
and fiux expected from a spherically symmetric stellar wind 
as a result of free-free emission. Combining optically thick 
and optically thin regimes they predict a radio-mm spec- 
trum with spectral index ~ -1-0.6. The fiux density expected 
from the stellar wind of the OB-type mass donor in Cyg X-1 
(assuming M ~ 2.5 x 10"*^ and Vini ~ 2000 km s"^) would 
be around 0.1 mjy at 100 GHz. Therefore we can see that 
neither the spectrum nor flux density are compatible with 
the 'standard' spherically symmetric stellar wind model. 

Another possibility is that the radio-mm spectrum is 
a combination of some emissive mechanism at radio wave- 
lengths, probably synchrotron, with a thermal component at 
(sub-)mm wavelengths. In the case that this thermal emis- 
sion arose in an optically thick dust cloud which peaked at a 
frequency of ~ 10^^ Hz (30/im), this corresponds to a tem- 
perature of ~ 150 K for the dust cloud. At a distance of 



2.5 kpc, a spherical cloud of radius > 3 x 10^^ cm would be 
required. This would easily enclose the entire binary system, 
and presumably significantly redden the colours of the OB 
companion star. In addition, such a large cloud would im- 
pose a minimum timescale for 50% variability of > 10 min. 
This cloud size is not unfeasible for a massive OB-type com- 
panion, although in order to be in thermal equilibrium at 
~ 150K the dust cloud would need to be much further from 
the star (at 10^"' cm from the star the equilibrium temper- 
ature is likely to stiU bo > lOOOK). 



4 CONCLUSIONS 

Cyg X-1 was previously known to have a fiat radio spec- 
trum from 2-15 GHz (Pooley et al. 1999) with a mean 
fiux density of ~ 14 mJy. This corresponds to an integrated 
synchrotron luminosity of 10^*^ erg s~^. A single previous 
observation at 250 GHz had implied that this fiat spec- 
trum extended to mm wavelengths (Altenhoff et al. 1994). In 
multiple simultaneous radio and mm observations we have 
confirmed the existence of a spectral component extending 
from cm through mm wavelengths with a very fiat spectrum 
and no evidence of either low- or high-frequency cut-offs. 
Furthermore, the likelihood that adiabatic expansion losses 
dominate in the emitting region shows that the generation of 
the outflow may be far more important to the energetics of 
accretion in Cyg X-1 than previously suspected. Presuming 
the emission to arise in a jet, and comparing luminosity and 
variability of this component with that from Cyg X-3 and 
GRS 1915-1-105, we infer that the outflow from Cyg X-1 is 
considerably steadier and has a significantly lower mass fiow 
rate. The radio-mm spectra of these X-ray binaries are much 
fiatter than those of the 'fiat-spectrum' AGN which gener- 
ally peak somewhere in the mm regime and fall off rapidly in 
the infrared. It is not at all clear whether models of partially 
self-absorbed synchrotron emission from conical jets which 
have been developed for these AGN can be extended to ap- 
ply to the much higher-frequency fiat-spectrum emission we 
observe from Cyg X-1, Cyg X-3 and GRS 1915+105. 

Detailed spectral measurements in the mm regime, 
preferably combined with simultaneous X-ray and radio ob- 
servations, should help us to improve the currently inad- 
equate understanding of the emission mechanisms in this 
unusual object. Measurement of the shortest variability 
timescale (expected to be < 1 sec for nonthermal emission 
with a brightness temperature of lO^^K, > 10 min for an 
optically thick dust cloud at 150 K, or > 1 hr for optically 
thin free-free emission) will be important in understanding 
the emissive mechanism. Equally important will be mea- 
surement of, or stringent upper limits to, the level of linear 
polarisation of the flat spectral component. Finally, we pre- 
dict that a flat spectrum (sub)mm component will also be 
detected from the persistent black hole candidate GX 339-4 
in the low/hard state. 
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